Best 

Available 

Copy 


AD-  774  472 


RESEARCH  DIRECTED  TOWARD  THE  USE 
OF  LONG  AND  INTERMEDIATE  PERIOD 
SEISMIC  WAVES  FOR  THE  IDENTIFICATION 
OF  SEISMIC  SOURCES 

Keith  McCamy 


Lamont  -  Doherty  Geological  Observatory 


Prepared  for: 


Advanced  Research  Projects  Agency 

Air  Force  Cambridge  Research  Laboratories 


October  1973 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


Lament -Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  Hew  York  1096*1 


CONTRACT  NO.  F19628-71-C-02*i5 


Project  No.  1795 

Task  No.  179500 

Work  Unit  No.  17950001 

SEMI-ANNUAL  TECHNICAL  REPORT  NO. 

October  1973 

CONTRACT  MONITOR:  Carlton  E.  Moll  neux 
Terrestrial  Sciences  Laboratory 


Approved  for  public  releast 


■  ■ 


4 


Sponsored  by 

Defense  Advanced  Research  Projects 
ARPA  Order  No.  i.795 
Monitored  by 

AIR  FORCE  CAMBRIDGE  RESEARCH  LAB 
R  FORCE  SYSTEMS  COMMA N 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS,  0 


-  ms  ■  -s&iSm-M 


’ibution  unlimit 


Unclassified 


Security  Cla»Kification  _____ 


DOCUMENT  CONTROL  DATA  .  R  &  D 

_ |5<Tiif<fy  cluu.lticatlon  o I  llllo,  body  cl  nhnimri  und  Indetlnd  imnotntlon  mu»l  be  anttrtd  when  the  overall  upon  Ih  cl.it/tUltd) 


i  ONlCtNA  TING  ACTIVITY  (Corporate  author)  20.  BE  PORT  IECUAITY  CLASSIFICATION 

Lamont-Doherty  Geological  Observatory  Unclassified 

Columbia  University  group  ’  “ 

Palisades,  New  York  10964 


i  mi  pom t  n  u  F. 

RELEARCh  DIRECTED  TOWARD  THE  USE  OF  LONG  AND  INTERMEDIATE  PERIOD 
SEISMIC  WAVES  I-'OR  TliE  IDENTIFICATION  OF  SEISMIC  SOURCES. 


*  DESCRIPTIVE  NOTES  (Type  ol  report  end  Inclusive  dotes) 


Scientific 


Interim 


5  authoriSI  (First  name,  middle  initial,  lost  name) 


Keith  McCamy 


P  REPORT  DATE 


7a.  TOTAL  NO.  OF  PACES 


lb.  NO  OF  REFS 


October  1973  10 


M.  CONTRACT  OR  GRANT  NO  OPtlCF  NO*  1(95  ORICNMOR'S  REPORT  NUMBER(S) 

F1962 8-71-0-024^  SEMI  ANNUAL  TECHNICAL  REPORT  NO.  4 

b.  PROJEC  T  NO  — • 

1795-00-01 
-  D0D  ELEMENT 
62710D 

d-  D0D  SUBELEMENT  n/a 


10  DISTRIBUTION  STATEMENT 


9b.  OTHER  REPORT  NO<S|  ( Any  other  numbers  that  moy  be  assigned 
this  report) 

AFCRL-TR-73-066 2 


A-  Approved  for  public  release;  distribution  unlimited 


II  SUPPLEMENTARY  notes 


112  SPONSORING  MILITARY  ACTIVITY 


This  research  was  supported  by  the  Air  Force  Cambridge  Research 
Defense  Advanced  Research  Project  Laboratories  (LW) 

Agency _ L.G.  Ilanscoin  Field 

'*  *BST"*CT  Bedford,  Maasachusettes  01730 _ 

Seismological  research  supported  by  Air  Force  Contract  F19628-71-C- 
024b  at  the  Lamont-Doherty  Geological  Observatory  is  summarized  for  the 
period  1  January  to  30  June  1973. 

"^During  this  period,  significant  advances  have  been  made  in  our 
understanding  of  the  physical  processes  of  natural  seismic  sources. 
Intra-plate  earthquakes  —  the  isolated  events  which  occur  in  the  rela¬ 
tively  seisinically  inactive  interiors  of  lithospheric  plates  —  have 
been  found  to  display  a  pattern  whose  consistency  on  a  global  scale  is 
emerging,  -‘■his  pattern  provides  important  constraints  on  the  mechanisms 
possible  for  driving  plate  tectonics  as  well  as  reducing  the  chances  of 
intra— plate  earthquakes  being  misidentified  as  explosive  sources.  The 
intra-plate  earthquake  swarm  at  Blue  Mt.  Lake,  N.Y.  provided  evidence 
that  helped  evoke  the  dilatancy  model  of  the  physical  process  of  rock 
failure  ^opening  the  door  to  the  studies  of  earthquake  sources,  usin^ 
str°ng-moUon  records  complement  theoretical  investigations  modelling 
the  dynamics  and  kinematics  of  shear  failure.  ° 

wavesThLrbeenSdefltaw^h1hi?h°!1H6eneitleS  °n  the  ProPaKation  of  surface 

obaervatlonally  by  obs^^ng  tSe  pass^fof^o^  "odelll"K 

nental  margin.  passage  of  Love  waves  across  a  conti- 
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MAJOR  SCIENTIFIC  ACCOMPLISHMENTS 


In  the  following  paragraphs,  scientific  accomplishments 
for  December  thru  June  1973,  are  summarised  following  the 
itemization  in  the  statement  of  work  of  this  contract. 


Line  -  Item  la 


Operation  of  the  Lamont-Doherty  network  of  three  long- 
period  and  intermediate-period  seismographic  stations  (the 
. al^sades ,  Sterling  Forest,  and  Ogdensburg)  has  continued 
during  the  last  six  months. 

Line  -  Item  lb 


Intra— plate  earthquakes,  are  fairly  common  in  North 
America  and  in  Asia.  It  is  these  lone  type  of  events  within 
areas  that  a -c  not  normally  thought  of  as  highl:  scismically 
active  that  must  be  identified  as  earthquakes  and  not  become 
a  false  alarm  by  misidentiflcation  as  a  possible  underground 
nuclear  explosion.  Many  of  these  events  appear  to  be  enriched 
in  the  high-frequency  part  of  the  spectrum.  A  preliminary 
analysis  indicates  that  they  fall  near  the  extremity  of  the 
earthquake  population  on  an  Ms  -  m  population  and  that  they 
do  not  fall  in  the  explosion  population. 

It  appears  likely  that  intra-plate  earthquakes  and  present- 
day  stresses  arc  closely  related  to  the  probable  drivinr 
mechanisms  of  plate  tectonics.  In  addition  to  being  of  para¬ 
mount  importance  to  geology,  understanding  such  mechanisms 
will  also  permit  global  generalizations  on  the  nature  of  intra- 

ra?oh?U?keS\  Syk°3  and  Sbar  havo  ctudled  the  focal  mechan- 
J"  ?  intra-plate  earthquakes  as  well  as  other  information 

on  in  situ  stress.  iheir  principle  result  is  that  a  large 
percentage  of  the  shocks,  particularly  those  not  located  near 
riid-ocean  ridges  or  behind  3ubduction  zones, are  characterized 
by  a  predominance  of  thrust  faulting.  In  the  few  areas  of  the 
worid  for  which  they  are  available,  other  determinations  of 
tne  state  of  stress  from  overcoring,  hydro fracturing,  and 
recent  crustal  movements  support  their  conclusion  that  larre 
areas  in  the  interiors  of  many  lithospheric  plat-s  arc  char¬ 
acterized  by  high  horizontal  compressive  stresses  that  appear 
to  be  uniform  in  orientation  over  sizable  regions  for  the  few 
cases  wnere  sufficiently  dense  data  arc  available.  Assuming 

a  t£ctonlr1n^iiP“t.hrn  °f  s,Te5n  ln  the  Plato  interiors  has 
origin ,  these  results  help  to  discriminate  between 
proposed  mechanisms  for  driving  the  plates.  For  example 
high  horizontal  compressive  stress  is  inconsistent  with  grav¬ 
itational  sinking  a3  the  dominant  motive  force  for  the  plates. 
n?SnSihl'?keS|ard  Sbar's  observations  are  compatible  with 
S  h  r?Ve?^°m  below  by  flow  in  thG  asthenosphere  or 

ridKC  croft-  aTheie  are  PUShin(!  from  tne  eeneral  area  of 
that  df  ?  ;  however,  several  focal  mechanisms 

that  do  not  fit  the  simplest  models  of  eravitatinnai  5lldln 


2. 


as  the  primary  motive  force  of  plate  tectonics.  Finally, 
the  spatial  distribution  of  the  directions  ol  principle 
horizontal  stress  shows  no  clear  relation  to  the  distribution 
of  hot  spots  which  argues  against  mantle  plumes  as  a  pro¬ 
posed  major  driving  mechanism.  Thus  their  study  favors  flow 
j.n  the  asthemsphere  as  the  most  likely  major  driving  mechanism 
of  plate  tectonics. 


Work  currently  in  progress  includes  developing  a  tectonic 
framework  for  the  observed  seismicity  in  central  Asia.  Both 
spatial  and  temporal  distribution  of  seismicity,  as  well  as 
fault  plane  solutions,  are  being  employed  toward  the  evolution 
of  a  large  scale  tectonic  model.  Previous  work  on  this  subject 
has  been  limited  to  larger  events,  generally  to  events  of  mh 
K^ater  than  5.5.  Recently,  techniques  have  been  developed0 
[Tat ham  and  Savino,  1973]  for  using  Very  Long  Period  Seismograph 
otations  to  discriminate  between  possible  focal  mechanisms  by 
observation  of  long-period  surface  waves.  This  technique  has 
been  applied  to  events  of  m  less  than  lf.0  at  distances  of 
n f °win,r'0  and  nas  been  expanded  to  use  the  long-period  components 
of  WM.oo  Stations.  Ane  method  is  presently  being  modified  for 
central.  Asia.  The  modification  involves  using  known 
fault  plane  solutions  to  smaller  events  in  adjoining  regions. 

h0r  of5erVhe  P°tential  for  greatly  expanding  data 
available  for  developing  a  tectonic  model  of  the  region,  as 

weli  as  offering  a  possible  tool  frr  discrimination  between 
observed  anomalous  events  and  underground  nuclear  explosions. 

Line  -  Item  lc 


_ 1:,  --UPture  zomt5  and  possible  locations  of  large  earthouakes 

plate . t  U  °  and  northern  boundary  of  the  Phllllplne  Sea - 

Mn,oL?I!:K3hal,l0W  earthquakes  (H  >  6.9)  that  have  occurred 
olnce  1?( d  have  been  relocated  to  determine  whenever  possible 

fupt^re  a^ng  the  northern  and  western  boundary 

e Jh  fP;  1  plnVea  plate.  Two  areas  along  the  seismic  zone 

qua  e°ii  ?n°  n0t  cxP^rlenced  a  large,  shallow,  earth¬ 

quake  in  at  least  30  years  and  are  areas  of  special  seismic 

potential.  A  compilation  of  iiistoric,  microscopic,  and  re- 

Juns  abo  ?Slt?n:nCT!;ar  d3t?  Shows  that  tho  PhiUlplne  fault  which 
runs  aboit  1200  km  from  Luzone  to  Mindinao  probably  has  ruptured 

over  at  least  half  its  length  during  the  last  105  years  Con- 
?lat|nnol  a?,t,ivlty  f1  the  PhlillPine  fault  Indicated  that 

accommodf  ted  by  t Ms  ?aSn?  *’e  Partlally  «  *ot.ll, 

.'he  Ryukyu  Island  arc,  between  123°  and  130°I.T  has  not 
experienced  a  large  shallow  earthquake  in  50  years ’and  is 
t®  ®  ve3y  aselsmic  for  a  plate  boundary.  \  study  of  historic 
Japanese  sources  indicates,  however,  that  large,  destructive 
earthquakes  with  associate  tsunamis  have  occurred  along  this 
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portion  of  the  plate  boundary  is  an  area  of  special  seismic 
potential. 

Line  -  Item  je 

Keceno  work  casting  much  light  on  the  actual  physical 
processes  taking  place  at  an  earthquake  hypocenter  has  been 
supported,  in  part,  by  this  contract.  Scholz,  Sykes,  and 
Aggarwal  have  developed  the  rock  dilatancy-water  diffusion 
model  of  the  earthquake  source  region  from  a  diverse  sei 
of  data,  some  of  which  was  obtained  from  the  seismograph 
network  around  the  Blue  Mt .  Like  earthquake  swarm.  This 
uilatarcy  model  goes  far  toward  explaining  the  nature  of  an 
earthquake  source  before,  during  and  after  failure.  They 
show  that  precursory  effects  such  as  crustal  movement,  changes 
in  seismic  velocities,  tilts,  fluid  pressure,  electrical  and 
magnetic  fields,  radon  emission,  and  frequency  and  distribution 
by  magnitude  of  small  local  shocks  occur  before  many,  perhaps 
all,  shallow  earthquakes  and  demonstrate  that  these  heretofore 
unrelated  premonitory  effects  oc:ur  at  a  characteristic  time 
before  the  earthquake  which  is  directly  relatablo  to  the  mag¬ 
nitude  of  the  event.  This  physical  basis  promises  to  lead 
us  soon  to  quantitative  deterministic  earthquake  prediction. 


Briefly,  the  model  describes  a  sequence  of  events  that 
occur  before  failure.  As  the  ambient  stress  rises  in  the 
source  region  due  to  tectonic  forces  the  rock  undergoes  an 
inelastic  increase  in  volume  produced  by  cracks  forming  and 
propagating.  This  phenomenon  is  called  dilatancy  and  can 
occur  at  stresses  as  low  as  half  the  breaking  strength  of  the 
rock.  While  cracks  are  forming  and  for  some  time  after,  the 
rock  will  be  undersaturated  because  of  the  finite  time  it 
tikes  water  to  diffuse  in  to  fill  the  new  crackr.  The  under- 


saouration  of  the  rocks  causes  the  propagation  velocity  of 
P  waves  to  decrease,  as  well  as  several  other  observed  pre¬ 
monitory  effects.  In  addition,  the  undersaturation  decreases 
the  pore  pressure  in  the  dilatant  region  and  hence  increases 
the  effective  stress,  a  phenomenon  called  dilatancy  hardening, 
which  actually  strengthens  the  rock,  inhibiting  rupture  despite 
a  continued  rise  in  ambient  stress.  As  fluid  flows  into  the 
dilatant  region  the  new  cracks  fill  and  the  rock  becomes  more 
saturated.  The  pore  pressure  cannot  begin  to  rise  until  the 
I’ock  is  wholly  saturated  but  when  it  does,  it  lowers  *-,->e 
effective  stress,  weakening  the  rock.  Since  the  tectonic 
stress  continues  to  rise  during  the  dilatant  period,  the 
rising  pore  pressure  triggers  the  earthouake  just  as  it  does 
through  fluid  injection  or  reservoir  filling.  Thus  dilatancy 
first  delays  the  earthquake  by  reducing  fluid  pressure  on 

t»Jthon  ^Peers  it  when  the  pore  pressure  is  recover- 
.  V  1  ^ime  scale  of  this  sequence  of  events  is  controlled 

Si.ih0rdJ£fU!!J?n  0f  fluid  and  consenuently  depends  on  the 
oise  of  the  dilatant  region,  relating  the  duration  of  the 
premonitory  stages  to  the  magnitude  of  the  earthquake. 


Tills  theory,  in  addition  to  promising  exciting  advances 
in  earthquake  prediciton,  provides  the  first  comprehensive 
physical  description  of  the  nature  of  the  earthquake  source 
re cion.  To  date,  the  dilatancy  model  seems  to  apply  best  to 
intraplate  events  where  the  dominant  stress  is  compmssional . 
extending  the  model  to  other  stress  regimes  is  an  active  area 
of  research. 

3.  Five  accelerograms  recorded  near  the  southeast  end 
of  the  rupture  of  June  27,  196 6  earthquake  at  Farkfield, 
^aliforn..a  have  been  doubly  integrated  to  give  ground  dis- 
placement  records.  Theoretical  ground  displacements  were 
caxculated  for  comparison  based  on  a  smoothly  propagating 
ramp  in  a  homogeneous  half  space.  Tiie  theoretical  disnlace- 
ment  agreed  reasonably  well  with  the  observations  of  horizon¬ 
tal  motion  before  the  arrival  of  the  surface  waves,  with  a 
best  fitting  propagation  velocity  of  2.8  to  3.0  km/sec.  A 
paper  describing  this  modelling  of  an  earthquake  source  is 
in  press  with  the  Bulletin  of  the  Seismological  Society  of 
wnenca  • 


Line  -  Item  If 


_  A  combined  theoretical  and  observational  study  of  the 

refraction  of  Love  waves  employs  a  method  developed  for  cal¬ 
culating  reflection  and  transmission  of  inhomogeneous  SII 
waves.  r  he  components  of  the  Love  waves  in  the  appronriate 

“etarnil?ed  to  obtain  the  transmission  and  reflection 
coefficients  for  Love  waves.  This  technique  agrees  excellently 

takesPr0'-1°^ofJ;nlte  el6?m,t  results.  A  typical  calculation 
takes  .0t.  seconds  on  an  IBM  3b0,  Model  91,  which  compares 

most  favorably  with  the  10  to  20  minutes  required  for  the 

and^n^t-f1"1?6  eiement  computation.  Love-wave  transmission 
and  refit ction  has  been  numerically  investigated  in  models 
a  continental  margin.  The  boundary  conditions  are  satis- 
^e?<eXJ\CtJy  at  the  lnterface  between  the  oceanic  and  the 
continental  model.  Diffracted  waves  for  the  problems  investi- 

5?  be  C?nS7dered  have  a  minor  influenceT  ^nJe  most 

the  energy  is  in  the  transmitted  or  reflected  Love  waves. 


Pre] 
amplitude 
OBS,  and 
Californi 
on  perioL 
the  trans 
observed 
to  the  oc 
longest  c 


iminary  observations  of  the  transmitted  Love  wave 
s  have  been  made  at  the  ocean  bottom  seismic  station, 
the  nearby  seismic  station  Byerly,  BKC,  in  Berkeley, 
a.  The  transmitted  amplitudes  are  strongly  dependent 
mi  t  tori  ProP®p^i°n  Prom  the  ocean  to  the  continent, 

amplitudes  are  largest  for  the  shortest  of  the 

Jw*  t°  SCC4*  For  Propagation  from  the  continent 
ean  the  transmitted  amplitudes  are  largest  for  the 
f  the  observed  period,  60  sec. 


To  explain  the  observed  trends  in  the  data  if  i  <=  goi-oH 
nan?l|aurder  the  ocean?  ^  ^  a  low"vel°oity  layer  in  the 

continental  Sn'SfSiM? 5  at,th?  »™««lng  of  the 
’  a°  bke  transmission  dependence  on 


angle  Incidence  Is  presently  under  investigation. 

?.  Richards  has  investigated  a  dynamic  model  of  earthquake 
faulting  to  help  q"antify  the  stress  drop  associated  with  an 
earthquake.  His  model  nucleates  at  the  fault  motion  at  a 
po4ot  and  allows  the  rupture  to  spread.  When  the  rupture  zone 
grows  over  a  large  part  of  the  fault  surface,  the  particle 
velocities  become  proportional  to  the  stress  drop.  The  pro¬ 
portionality  constant  may  vary  for  a  variety  of  earth  material. 
Westerly  granite,  for  example,  gives  60  cm/sc c  particle 
velocity  for  each  100  bars  of  stress  drop.  A  paper  describing 
this  work  3s  in  press  with  the  Journal  of  Applied  Physics. 

3.  A  major  problem  in  the  development  of  a  theory  for 
seismic  body  waves  lies  in  the  fact  that  P  and  S  waves  do 
not  propagate  independently  in  a  heterogeneous  medium  like 
the  earth.  Richards  has  approached  this  problem  in  terms  of 
potentials  for  P  and  S  v/aves,  and  it  has  been  found  possible 
to  display  a  coupling  coefficient  between  P  and  S7  motion. 

The  method  permits  many  known  solutions  for  scalar  waves  in 
the  heterogeneous  media  to  be  adapted  to  the  study  of  P  waves 
in  seisnology.  A  manuscript  describing  this  potential 
method  is  nearing  completion,  with  submission  to  the  Bull, 
oeis.  Soc .  of  Amer.  e/nected  before  the  end  of  September. 
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